


® BUIMDESREPUBLIK © Off nl gungsschrift 
DEUTSCHLAND DE 10136754A1 




DEUTSCHES 
PATENT- UND 
MARKENAMT 



® Aktenzeichen: 101 36 754.6 

@ Anmeldetag: 27. 7.2001 

® Offenlegungstag: 13. 2.2003 



® Int. CI. 7 : 

G 01 N 9/24 

G 01 F 23/296 
G01 N 29/16 
G 21 C 17/022 



CO 



UJ 

Q 



@ Anmelder: 

Siemens AG, 80333 Munchen, DE 



® Erfinder: 

Bechtold, Mario, 91334 Hemhofen, DE; Vester, 
Markus, Dr., 90471 Nurnberg, DE 



J6) Entgegenhaltungen: 


DE 


35 41613 1 


US 


55 54 936 


US 


52 49 463 


US 


47 86 857 


US 


3812 422 



##0Bechtold, Mario, 91334 Hemhofen, DE; Vester, 
Markus, Dr., 90471 Nurnberg, DE 

H 



15 

L 



Die folgenden Angaben sind den vom Anmelder eingereichten Unterlagen entnommen 

Prufungsantrag gem. § 44 PatG ist gestellt 

® Verfahren und Vorrichtung zur Dichtebestimmung 

® Bei dem Verfahren und der Vorrichtung zur Dichtebe- 
stimmung eines in einem Behalter (20) befmdlichen Me- 
diums (30) mittels einer in den Behalter (20) hineinragen- 
den Sonde (10) wird ein Abfragesignal (SI ) in die Sonde 
(10) eingespeist und ein durch das Abfragesignal (S1) her- 
vorgerufenes Antwortsignal (S2) detektiert. Mittels des 
Zeitverlaufs des Antwortsignals (S2) wird ein Reflexions- 
faktorverlauf (R) des Mediums (30) langs der Sonde (10) 
bestimmt. Der Reflexionsfaktorverlauf (R) wird zur Ermitt- 
lung des Dichteverlaufs (p) des Mediums (30) langs der 
Sonde (10) verwendet. 
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Beschreibung 

[0001] Die Erfindung bezieht sich auf ein Verfahren und eine Vorrichtung zur Dichtebestimmung eines in einem Be- 
halter befindlichen Mediums mittels einer in den Behalter hineinragenden Sonde. 
5 [0002] Wahrend des Bctriebes eines Kernkraftwerks oder auch wahrend verschiedener Prozessablaufe in der chemi- 
schen Verfahrenstechnik ist es notwendig, eine MessgroBe innerhalb eines befiillten Behalters zu ermitteln. Aufgrund der 
hohen Temperatur oder auch der gegebenen falls innerhalb des Behalters herrschenden radioaktiven Strahlungsintensitat 
konnen innerhalb des zu iiberwachenden Behalters keine elektronischen Messmittel eingesetzt werden. Insbesondere 
beim Betrieb eines Kernreaktordruckbehalters ist die Kenntnis uber das Vbrhandensein und den Zustand des im Reaktor- 
10 behalter befindlichen Kiihlmediums von besonderem Interesse. Diese Messinformation soil auch in einem Storfall zuver- 
lassig ermittelt werden konnen. 

[0003] Aus der DE 35 4 1 6 1 3 Al ist ein Verfahren zur Fullstandsmessung in einem Reaktordruckbehalter bekannt Bei 
diesem Verfahren wird ein DifFerenzdruck zwischen dem eigentlichen Reaktordruckbehalter und einer Vergleichssaule 
gemessen. Daraus und unter Beriicksichtigung der Flussigkeits- und Dampfdichten wird in einer Auswerteeinheit der 
15 Fullstand ermittelt. Dieses Messverfahren weist jedoch insbesondere bei den in einem Storfall vermehrt auftretenden 
Turbulenzen innerhalb des Kiihlmediums im Reaktordruckbehalter Ungenauigkeiten auf. 

[0004] In dem Fachaufsatz von H. MuUer-Nordhoff "Fullstandsmessung nach der TDR-Methode im Aufbereitungs- 
und kohlenveredelungsbetrieben" aus Gluckauf-Forschungshefle 44 (1983), Heft 5, Seiten 214 bis 218 ist ein anderes 
Messverfahren zur Fullstandsmessung beschrieben. Dieses sogenannte Zeitbereichs-Reflektometrie- Verfahren (Time 
20 Domain Reflectometry = TDR) beruht auf der Messung der Laufzeit eines Hochfrequenzimpulses langs einer Sonde bis 
hin zu einer Unstetigkeitsstelle, wo ein Teil der eingespeisten Hochfrequenzenergie reflektiert wird. Diese Unstetigkeits- 
stelle wird dann als Fullstandshohe des im Behalter befindlichen Mediums interpretiert. 

[0005] Dieses auf einer reinen Laufzeit-Messung basierende TDR- Verfahren zur Fullstandsmessung wurde in seiner 
Messgenauigkeit dadurch verbessert, dass an der in das zu detektierende Medium hineinragenden Sonde an definierter 
25 Stelle Referenzmarken angebracht wurden. Dadurch erhalt man im Antwortsignal zusatzliche durch die Referenzmarken 
bedingte Reflexionsimpulse, deren Reflexionsorte im Gegensatz zu der zu detektierenden Fullstandshohe exakt bekannt 
sind. Deshalb konnen die fur die Referenzmarken detektierten Laufzeiten zur Kalibrierung des TDR- Verfahrens heran- 
gezogen werden. Ein solches Verfahren ist in der US 5,249,463 beschrieben. 

[0006] Die auf einer Laufzeitmessung basierenden Standard-TDR- Verfahren weisen aber bei einem Einsatz zur FU11- 

30 standsmessung eines fliissigen Kiihlmediums mit zunehmender Temperatur eine steigende Messungenauigkeit auf. Mit 
steigender Temperatur verdampft namlich zum einen ein Anteil des Kiihlmediums und zum anderen nahem sich die Aus- 
breitungsgeschwindigkeiten in dem dampffbrmigen und dem fliissigen Anteil des Mediums einander an, so dass keine 
exakt definierte Unstetigkeitsstelle beispielsweise in Form einer Flussigkeitsoberflache mehr vorliegt. Bei Wasser steigt 
insbesondere oberhalb von 200°C die Dichte des Wasserdampfs an und die des fliissigen Wassers fallt, bis sich beide 

35 Werte im kritischen Punkt bei 374°C treffen. Mangels markanter Unterschiede in der Ausbreitungsgeschwindigkeit fuh- 
ren die Standard-TDR- Verfahren in dem genannten hohen Temperaturbereich also zu einem nicht unerheblichen Mess- 
fehler, der zusatzlich durch TVirbulenzen und die damit verbundene Bildung von Mischphasen aus dem dampfiformigen 
und dem fliissigen Aggregatszustand des Mediums erhoht wird. Der Aggregatszustand wird manchmal auch als Phasen- 
zustand oder einfach als Phase bezeichnet 

40 [0007] Zur Umgehung dieser Schwierigkeiten wird in der US 4,786,857 ein modifiziertes TDR- Verfahren beschrieben, 
bei dem anstelle einer einzigen langen Sonde eine Vielzahl deutlich kiirzerer Messsonden innerhalb des zu iiberwachen- 
den Behalters vorgesehen sind. Geeignete Ausfuhrungsformen dieser kurzen Messsonden sind in der US 5,554,936 be- 
schrieben, Jede dieser etwa 30 cm langen Einzelsonden ist mit einer eigenen TDR-Auswerteeinheit verbunden. Aus den 
Reflexionsstellen am Anfang und am Ende jeder Einzelsonde lasst sich mit den Standard-TDR-Methoden eine mittlere 

45 Ausbreitungsgeschwindigkeit im Messbereich jeder Einzelsonde ermitteln. Zusatzlich werden vorab ermittelte Antwort- 
signale fur verschiedene Betriebszustande der Messsonde (Messsonde zum Teil in fliissigem Wasser eingetaucht, Mess- 
sonde in Wasserdampf, Messsonde in Mischphase aus Wasserdampf und fliissigem Wasser) bei der Auswertung mit be- 
riicksichtigt. Neben einer Erhohung der Messgenauigkeit lassen sich mit diesem Verfahren auch zusatzliche Informatio- 
nen uber das Medium in dem Behalter ermitteln. Allerdings ist der apparative Aufwand aufgrund der vielen Einzelson- 

50 den und ihrer zugehorigen Auswerteeinheiten betrachtlich. Weiterhin wird die Ortsauflosung und damit die Messgenau- 
igkeit dieses modifizierten TDR- Verfahrens durch die Anzahl und die Lange der Einzelsonden begrenzt. 
[0008] Die Aufgabe der Erfindung besteht deshalb darin, ein Verfahren und eine Vorrichtung zur Dichtebestimmung 
eines in einem Behalter befindlichen Mediums mittels einer in den Behalter hineinragenden Sonde anzugeben, welche 
eine hohe Messgenauigkeit bei gleichzeitig geringem apparativen Aufwand ermoglichen. 

55 [0009] Zur Losung dieser Aufgabe wird ein Verfahren entsprechend den Merkmalen des unabhangigen Patentan- 
spruchs 1 sowie eine Vorrichtung entsprechend den Merkmalen des unabhangigen Patentanspruchs 13 angegeben. 
[0010] Bei dem erfindungsgemafien Verfahren zur Dichtebestimmung eines in einem Behalter befindlichen Mediums 
mittels einer in den Behalter hineinragenden Sonde wird 

60 a) ein Abfragesignal in die Sonde eingespeist, 

b) ein durch das Abfragesignal hervorgerufenes Antwortsignal detektiert, 

c) ein Reflexionsfaktor- oder Impedanzverlauf des Mediums langs der Sonde mittels des Zeitverlaufs des Antwort- 
signals und 

d) ein Dichteverlauf des Mediums langs der Sonde mittels des Reflexionsfaktor- oder Impedanzverlaufs bestimmt. 

65 

[0011] Die erfindungsgemSBe Vorrichtung zur Dichtebestimmung eines in einem Behalter befindlichen Mediums mit- 
tels einer in den Behalter hineinragenden Sonde umfasst mindestens 
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a) Mittel zur Erzeugung und Einspeisung eines Abfragesignals in die Sonde, 

b) Mittel zur Detektion und Auswertung eines durch das Abfragesignal hervorgerufenen Antwortsignals, wobei 

c) die Mittel zur Detektion und Auswertung ausgebildet sind zur Bestimmung eines Reflexionsfaktor- oder Impe- 
danzverlaufs des Mediums langs der Sonde mittels des Zeitverlaufs des Antwortsignals und zur Bestimmung eines 
Dichteverlaufs des Mediums langs der Sonde mittels des Reflexionsfaktor- oder Impedanzverlaufs. 5 

[0012] Die Erfindung beruht dabei auf der Erkenntnis, dass eine hohe Messgenauigkeit auch dann zu erreichen ist, 
wenn nur eine einzige Messsonde wie bei einem Standard-TDR-Verfahren zur Bewachung des in dem Behalter befind- 
lichen Mediums verwendet wird. Im Gegensatz zu einem Standard-TDR-Verfahren wird jedoch nicht eine einfache 
Laufzeitmessung eines an einer Unstetigkeitsstelle reflektierten Antwortsignals erfasst und ausgewertet, sondem das to 
Antwortsignal mit seinem kompletten zeitlichen Verlauf fiir die Auswertung herangezogen. Der Informationsgehalt des 
Antwortsignals wird somit deutlich besser ausgenutzt als bei einer blofien Auswertung einzelner Reflexionsimpulse in- 
nerhalb des Antwortsignals. Es wurde erkannt, dass unter Berucksichtigung auch der anderen Tfeile des Zeitverlaufs des 
Antwortsignals eine ortliche Verteilung der Impedanz in Richtung der Sonde und daraus einen Dichteverlauf des Medi- 
ums langs dieser Richtung ermittelt werden konnen. 15 
[0013] Die Impedanz, die das eingespeiste Abfragesignal bei seinem Weg durch die Sonde erfahrt, wird zum einen 
durch die vorgegebene Sonden-Geometrie und zum anderen durch das die Sonde umgebende Medium bestimmt. Fiir die 
Auswertung ist der zuletzt genannte Anteil maBgeblich. Alternativ zur Bestimmung des Impedanzverlaufs langs der 
Sonde kann auch der Reflexionsfaktorverlauf in Sondenrichtung ermittelt werden. Die elektrische Impedanz und der Re- 
flexionsfaktor lassen sich iiber eine wohlbekannte einfache Formelbeziehung ieicht ineinander uberftihren. Beide Dar- 20 
stellungsformen enthalten den gleichen Informationsgehalt uber das zu uberwachende Medium und sind deshalb aquiva- 
lent einsetzbar. 

[0014] Das Medium kann in Form eines Stoffes, eines Stoffgemisches oder auch einer Schichtung mehrerer Stoffe 
(beispielsweise Wasser am Behalterboden und dariiber normale Umgebungsluft) ausgebildet sein. AuBerdem ist es auch 
moglich, dass das zu uberwachende Medium gleichzeitig in mehreren Aggregatszustanden vorliegt. Die einzelnen Ag- 25 
gregatszustande (beispielsweise fliissig und gasfbrmig) konnen dabei langs der Sonde in beliebig abwechselnder Reihen- 
folge vorkommen. Dariiber hinaus ist auch eine Mischphase aus beiden Aggregatszustanden moglich. Obwohl sich ge- 
rade im zuletzt genannten Fall keine ausgepragten Unstetigkeitsstellen im Reflexionsfaktor- oder Impedanzverlauf langs 
der Sonde ausbilden und somit ein Standard-TDR-Verfahren nicht mehr anwendbar ist, sind dennoch Informationen uber 
den Zustand des Mediums aus der brtlichen Verteilung dieser wichtigen KenngrdBen langs der Sonde zu entnehmen. 30 
[0015] Insbesondere lasst sich auch der Dichteverlauf aus dem Verlauf der elektrischen Impedanz und des Reflexions- 
faktors ermitteln. Die Dichte ist eine wichtige KenngroBe zur Beurteilung der aktuellen Verhaltnisse innerhalb des zu 
iiberwachenden Behalters. 

[0016] Die Messgenauigkeit kann bei dem Verf ahren auf einfache Weise iiber die Anzahl der berucksichtigten Abtast- 
werte des Abfragesignals erhoht werden. Eine Erhohung des apparativen Aufwandes im Innenbereich des Behalters in 35 
Form von zusatzlichen Messsonden wie beim Stand der Technik ist also nicht erforderlich. Damit lasst sich das Verf ah- 
ren auch an einer bereits bestehenden tJberwachungseinrichtung mit einer einzigen im Inneren des Behalters angebrach- 
ten Sonde nachrusten, wobei zugleich eine hohere Messgenauigkeit erreicht wird. Die wesentlichen Anderungen hierfur 
konnen in der leichter zuganglichen Auswerteeinheit vorgenommen werden. 

[0017] Vorteilhafte Ausfiihrungsformen des Verfahrens und der Vorrichtung ergeben sich aus den von Anspruch 1 bzw. 40 
Anspruch 13 abhangigen Anspruchen. 

[0018] Bei einer ersten vorteilhaften Ausfuhrungsvariante wird der Reflexionsfaktor- oder Impedanzverlauf mittels ei- 
nes Rekursionsalgorithmus aus dem Zeitverlauf des Antwortsignals bestimmt. Dadurch konnen auch gegebenenfalls auf- 
tretende Mehrfachreflexionen zwischen Teilbereichen der Messsonde mit erfasst und aus dem Antwortsignal herausge- 
rechnet werden. w 45 

[0019] Bevorzugt ist ein Verfahren zur Ermittlung des Reflexionsfaktor- und Impedanzverlaufs, bei dem in der Aus- 
werteeinheit eine fiktive Einteilung des Mediums in mehrere Schichten, die quer zur Sonde verlaufen, erfolgt. Die zu- 
nachst noch unbekannte raumliche Ausdehnung in Langsrichtung der Sonde ist jeweils so bemessen, dass die Laufzeit 
des Abfragesignals in jeder Schicht gleich ist. Beginnend mit der in Einspeiserichtung des Abfragesignals als erstes an- 
geordneten Schicht wird fur jede Schicht ein aktueller Reflexionsfaktor, eine aktuelle Ausbreitungsgeschwindigkeit so- 50 
wie die aktuelle Schichtausdehnung in Einspeiserichtung des Abfragesignals ermittelt. Durch Zuordnung des aktuellen 
Reflexionsfaktors zu dem aktuellen Ort ergibt sich sukzessive die gewiinschte ortliche Verteilung des Reflexionsfaktor- 
verlaufs. Bei Bedarf kann diese Verteilung auch in eine Impedanzverteilung umgerechnet werden. Der aktuelle Reflexi- 
onsfaktor wird dabei der Grenze zwischen der gerade betrachteten aktuellen Schicht und der angrenzenden im vorherge- 
henden Rekursionsschritt des betrachteten Schicht zugeordnet. Die Ermittlung erfolgt unter Berucksichtigung des emp- 55 
fangenen und gespeicherten zeitlichen Verlaufs des Abfragesignals, des bekannten eingespeisten Antwortsignals und der 
in den vorherigen Schritten des Rekursionsalgorithmus ermittelten Reflexionsfaktoren an den Grenzen zwischen den 
vorhergehenden Schichten. Die aktuelle Ausbreitungsgeschwindigkeit der aktuellen Schicht wird aus dem aktuellen Re- 
flexionsfaktor sowie der in dem vorhergehenden Rekursionsschritt berechneten Ausbreitungsgeschwindigkeit der an- 
grenzenden vorhergehenden Schicht bestimmt. " go 
[0020] Die Ausbreitungsgeschwindigkeit kann auch als Mittelwert der sich innerhalb der aktuellen Schicht einstellen- 
den Geschwindigkeitsverteilung interpretiert werden. Da die Schichtabmessungen jedoch iiblicherweise sehr klein ge- 
wahlt werden konnen, stellt dies keine Beeintrachtigung der Messgenauigkeit des Verfahrens dar. 
[0021] Mit Kenntnis der aktuellen Ausbreitungsgeschwindigkeit kann der innerhalb der fest vorgegebenen Laufzeit 
zuriickgelegte Weg und damit die raumliche Ausdehnung der aktuellen Schicht ermittelt werden. Dem aktuellen Refle- 65 
xionsfaktor wird der Ort der Grenze zwischen aktueller und angrenzender, im vorherigen Rekursionsschritt betrachteter 
Schicht zugewiesen. Dieser Grenzort ist im vorhergehenden Rekursionsschritt uber die ermittelte Ausdehnung der vor- 
hergehenden Schicht errechnet worden. Nach Durchlaufen aller Rekursionsschritte erhalt man die Ortliche \ferteilung des 
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Reflexionsfaktors langs der Sonde, 

[0022] Die Zahl der Schichten, in die der Messbereich langs der Sonde eingeteilt wird, ist gleich der um 1 reduzierten 
Anzahl der aquidistanten Zeitintervalle, in die das Antwortsignal unterteilt wind. Die sich im Laufe des Rekursionsalgo- 
rithmus ergebenden lokalen Abtastwerte des Reflexionsfaktor- oder Impedanzverlaufs sind dagegen in der Regel nicht 
aquidistant voneinander beabstandet. Dennoch lasst sich auch die ortliche Auflosung einfach durch Erhohung der zeitli- 
chen Abtastrate verbessem. 

[0023] Vorzugsweise wird zur Dichteermittlung zunachst in jedem Rekursionsschritt eine aktuelle relative Dielektrizi- 
tatskonstante fur die gerade betrachtete aktuelle Schicht bestimmt. Dies geschieht unter Zugrundelegung des fiir die ak- 
tuelle Schicht bestimmten aktuellen Reflexionsfaktors und der ftir die im vorherigen Rekursionsschritt bestimmten rela- 
tiven Dielektrizitatskonstante der angrenzenden vorhergehenden Schicht. Aus der aktuellen relativen Dielektrizitatskon- 
stante wird ein aktueller Dichtewert fur die gerade betrachtete Schicht errechnet. Insbesondere bei einer vorgegebenen 
oder angenomxnenen Temperatur im Behalterinneren besteht namlich ein eindeutiger Zusammenhang zwischen relativer 
Dielektrizitatskonstante und der Dichte des im Behalter befindlichen Mediums. 

[0024] Bevorzugt ist eine weitere Variante, bei der aus der Dichte eine weitere interessante KenngroBe wie z. B. ein 
Fullstand innerhalb des Behalters oder auch ein Aggregatszustand ermittelt wird. Bei Zugrundelegung entsprechend 
kombinierter Formelbeziehungen lassen sich diese zuletzt genannten KenngroBen jedoch auch direkt aus dem Reflexi- 
onsfaktor- oder Impedanzverlauf bestimmen. 

[0025] Bevorzugt isl eine weitere Ausfuhrungsform, bei der zusatzlich eine Temperatur des Behalterinneren detektiert 
und bei der Auswertung mit beriicksichtigt wird. Je genauer die Temperatur innerhalb des Behalters bekannt ist, um so 
exakter lasst sich der Dichteverlauf langs der Sonde aus dem Reflexionsfaktor- oder Impedanzverlauf bestimmen. Des- 
halb ist es gunstig, wenn langs der Messsonde mehrere Temperaturmessstellen vorgesehen sind. Dann kann fur die Be- 
rechnungen in den jeweiligen Schichten ein moglichst genauer Wert der lokalen Temperatur verwendet werden. FUr eine 
ausreichende Messgenauigkeit ist es jedoch nicht unbedingt erforderlich, innerhalb jeder Schicht eine gesonderte Tem- 
peraturmessung durchzufuhren. Oblicherweise treten namlich keine signifikanten Temperaturschwankungen im Bereich 
der typischerweise verwendeten Schichtabmessungen.(z. B. einige cm) auf. 

[0026] Die Messgenauigkeit wird auBerdem im Rahmen einer weiteren Variante dadurch gesteigert, dass eine oder 
mehrere Referenzmarken an definierter Stelle der Sonde vorgesehen sind. Diese Referenzmarken fiihren zu besonders 
markanten Reflexionsimpulsen im Abfragesignal und bilden sich auch in dem errechneten Reflexionsfaktor- oder Impe- 
danzverlauf ab. Da die Position der Referenzmarken bekannt ist, kann somit eine Korrektur oder Kalibrierung des Ver- 
fahrens vorgenommen werden. Eine weitere Uberpriifungsmoglichkeit ergibt sich aus der Reflexion am Ende der Sonde. 
Aus der Laufzeit dieser Endreflexion und weiteren bekannten Vorgaben, beispielsweise der anfanglich vorhandenen 
Fiillmenge des Mediums, kann eine mittlere Ausbreitungsgeschwindigkeit langs der kompletten Sonde bestimmt wer- 
den. Damit lassen sich die in den einzelnen Schichten errechneten Werte auf Plausibilitat uberprufen und gegebenenfalls 
entsprechend anpassen. Neben der Laufzeit der Endreflexion kann auch deren Amplitude ausgewertet werden. Die Am- 
plitude gibt iiber die Absorption und die Leitfahigkeit im Medium Aufschluss. Daraus lasst sich beispielsweise ableiten, 
ob die Ionenkonzentration im Medium gestiegen ist. Ein solcher Anstieg ware ein Indiz auf das \brhandensein uner- 
wiinschter dissoziierter Zusatze, z. B. von Borsaure, oder auf die Einwirkung einer ionisierenden Strahlung, z. B. einer 
radioaktiven Strahlung. 

[0027] Bevorzugt hat das Abfragesignal die Form eines sprungformigen Signals, eines impulsformigen Signals oder 
eines Chirp-Signals, Andere Signalformen sind jedoch grundsatzlich ebenfalls moglich. Gunstigerweise hat das Abfra- 
gesignal dariiber hinaus eine Bandbreite von mindestens 50 MHz, vorzugsweise von mindestens 100 MHz. Die verwen- 
dete Bandbreite bestimmt neben der zeitlichen Abtastrate die Ortsauflosung und damit die Messgenauigkeit des \ferfah- 
rens. So ergibt sich bei einem Abfragesignal mit einer typischerweise eingesetzten Bandbreite von 1 GHz eine Ortsauf- 
losung in der GroBenordnung von 7,5 cm. Bei einer sprunghaften Anderung des Dichteverlaufs beispielsweise an einer 
Grenzflache des Mediums kann jedoch durch Berucksichtigung von ZusatzmaBnahmen wie etwa einer Kurvenanpassung 
("Fit) einer bekannten Sprungantwort an die gemessenen Daten oder eines Verfahrens zur hochauflosenden Entfaltung 
beispielsweise mit einem Wienerfilter eine deutlich hohere Messauflosung in der GroBenordnung von unter 1 cm erreicht 
werden. Die Messauflosung kann in diesem Fall auch durch ein Extrapolationsverfahren im Frequenzbereich oder eine 
Zusatzmessung im Zeitbereich zur genauen Ermittlung der zeitlichen Lage des Impedanz- oder Dichtesprungs erhoht 
werden. 

[0028] Bevorzugte Ausfuhrungsbeispiele werden nunmehr anhand der Zeichnung naher erlautert, Zur Verdeutlichung 
ist die Zeichnung nicht maBstablich ausgefuhrt, und gewisse Merkmale sind schematisiert dargestellt. Im Einzelnen zei- 
gen 

[0029] Fig. 1 eine Vorrichtung zur ortsaufgelosten Reflektionsfaktor- und Dichtebestimmung eines Mediums, 
[0030] Fig. 2 eine Schichteinteilung des Untersuchungsgebiets der Vorrichtung nach Fig. 1 sowie orts- und zeitaufge- 
loste Darstellung der das Antwortsignal bildenden Teilsignale, 

[0031] Fig. 3 einen Zusammenhang zwischen relativer Dielektrizitatskonstante und Dichte bei Wasser, 
[0032] Fig. 4 eine in der Vorrichtung gemaB Fig. 1 eingesetzte Sonde, 
[0033] Fig. 5 eine Sonde mit Temperaturmessstellen und 
[0034] Fig. 6 eine Sonde mit Referenzmarken. 

[0035] Einander entsprechende Teile sind in den Fig. 1 bis 6 mit denselben Bezugseichen versehen. 
[0036] In Fig. list eine Vorrichtung zur ortsaufgelosten Reflektionsfaktor- und Dichtebestimmung eines Mediums 30, 
das sich in einem Behalter 20 befindet, dargestellt. Ein in einer Generatoreinheit 4 erzeugtes elektrisches Abfragesignal 
SI wird iiber ein Koppelelement 5 und eine Vorlaufstrecke 15 in eine Sonde 10 eingespeist, die in dem Behalter 20 an- 
geordnet ist. Das zu detektierende Medium 30 ist im dargestellten Beispielsfall zweischichtig. Im unteren Teil des Be- 
halters 20 befindet sich eine WasserbefiJUung 31, im oberen Teil des Behalters 20 eine Luft- oder Wasserdampfbefiillung 
32. Das Abfragesignal S 1 wird in Abhangigkeit des die Sonde 10 umgebenden Mediums 30 reflektiert und lauft als Ant- 
wortsignal S2 iiber die Vorlaufstrecke 15 und das Koppelelement 5 zu einer Auswerteeinheit 3, die zusammen mit der 
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Generatoreinheit 4 in einer Sende-ZEmpfangseinheit 2 angeordnet ist, zuriick. In der Auswerteeinheit 3 wird das Ant- 
wortsignal S2 detektiert und aus seinem zeitlichen Verlauf eine ortliche Verteilung eines Reflexionsfaktors R fur das Me- 
dium 30 in Langsrichtung der Sonde 10 ermittelt. Diese Langsrichtung ist in Fig. 1 durch die Ortskoordinate d angezeigt. 
[0037] Anhand des ermittelten Reflexionsfaktorverlaufs R konnen in der Auswerteeinheit 3 weitere fur eine Zustands- 
bestimmung des Mediums 30 interessante KcnngroBcn abgeleitet werden. Anstelle des Verlaufs des Reflexionsfaktors R 5 
kann beispielsweise auch ein Verlauf einer elektrischen Impedanz Z des Mediums 30 langs der Sonde 10 angegeben wer- 
den. Der Reflexionsfaktor R und die Impedanz Z lassen sich ohne weiteres ineinander umrechnen. 
[0038] Aus dem Reflexionsfaktorverlauf R oder dem aquivalenten Impedanzverlauf Z wird als weitere MessgroBe M 
der Verlauf der Dichte p des Mediums 30 ermittelt. Weiterhin lassen sich Informationen iiber den Fullstand oder den Ag- 
gregatszustand des Mediums 30 bestimmen. l0 
[0039] Im Gegensatz zu den bekannten Standard-TOR(Time Domain Reflectometry)- Verfahren wird bei dem der Vor- 
richtung nach Fig. 1 zugrunde liegenden Verfahren nicht die Laufzeit des Abfragesignals S 1 bis zu einer bestimmten Re- 
flexionsstelle - im Beispiel von Fig. 1 beispielsweise der Grenzflache zwischen der WasserbefuUung 31 und der Luft- 
oder Dampfbefullung 32 - ermittelt. Vielmehr wird der komplette Verlauf des Reflexionsfaktors R oder der Impedanz I 
langs der Sonde 10 ermittelt. Damit steht in der Auswerteeinheit 3 ein deutlich hoherer Informationsgehalt zur VerfU- 15 
gung. Insbesondere ist eine MessgroBenermittlung auch dann moglich, wenn im Behalter 10 - anders als im in Fig. 1 dar- 
gestellten Fall - keine scharf definierte Reflexionsstelle in Form einer Grenzflache zwischen unterschiedlichen Stoffen 
vorhanden ist. Gerade dann kommt es bei Standard-TDR- Verfahren zu Messungenauigkeiten. Eine derartige Situation 
kann sich aber durchaus einstellen, insbesondere wenn es sich bei dem Behalter 10 um einen Kernreaktordruckbehalter 
handelt, der zur Kuhlung mit Wasser 31 befullt ist. Im Storfall konnen sich aufgrund der dann auftretenden hohen Tern- 20 
peratur und der TWbulenzen im Medium 30 ohne weiteres Mischphasen zwischen dem Wasser 31 und dem Wasserdampf 
32 ausbilden. Bei hohen Temperaturen nahem sich die elektrische Impedanz und damit auch die Dichte von flussigern 
Wasser und von Wasserdampf einander an. Im Storfall ist es deshalb durchaus moglich, dass keine mit einem Standard- 
TDR- Verfahren uber eine Laufzeitmessung langs der Sonde 10 detektierbare scharfe Reflexionsstelle mehr vorliegt. 
[0040] Das Verfahren zur Bestimmung des Reflexionsfaktorverlaufs R ist dagegen auch in dieser Grenzsituation funk- 25 
tionstiichtig und ermdglicht eine Aussage iiber den Zustand im Behalterinneren. Zur Ermittlung des Reflexionsfaktorver- 
laufs R wird in der Auswerteeinheit 3 ein eindimensionaler Rekursionsalgorithmus unter Beriicksichtigung des Abfrage- 
signals SI und des Antwortsignals S2 durchgefuhrt. Das Funktionsprinzip dieses Rekursionsalgorithmus wird im Fol- 
genden anhand der Darsteilungen von Fig. 2 erlautert. 

[0041] Als Abfragesignal S 1 wird eine Sprungfunktion in die Sonde 10 eingespeist. Ihr Verlauf uber der Zeit t sowie 30 
ein beispielhafter zeitlicher Verlauf des empfangenen Antwortsignals S2 ist in dem oberen bzw. dem unteren Diagramm 
von Fig. 2 gezeigt. Zur Auswertung werden die Sonde 10 und das sie umgebende Medium 30 zunachst in mehrere fiktive 
Schichten 101, 102 und 103 eingeteilt. Die Einteilung wird dabei so vorgenommen, dass das Abfragesignal SI fur das 
Passieren (Hin- und Ruckrichtung) jeder der Schichten 101, 102 und 103 eine identische fiktive Laufzeit X benotigt. Da 
die Ausbreitungsgeschwindigkeiten innerhalb der einzelnen Schichten 101, 102 und 103 in der Regel nicht identisch 35 
sind, ergeben sich bei Vorgabe jeweils gleicher Laufzeiten T unterschiedliche raumliche Ausdehnungen der einzelnen 
Schichten 101, 102 und 103 in Richtung der Ortskoordinate d ein. Diese raumlichen Ausdehnungen sind zunachst unbe- 
kannt. 

[0042] Das Abfragesignal SI und das Antwortsignal S2 werden zeitlich aquidistant mit einem RastermaB, das genau 
der Laufzeit T entspricht, abgetastet. Sowohl die Schichteinteilung der Sonde 10 als auch die zeitliche Rasterung sind in 40 
Fig. 2 kenntlich gemacht. Der zeitlichen Abtastung ist die diskrete Zeitkoordinate k zugeordnet, der ortlichen Schicht- 
einteilung der Sonde 10 die diskrete Ortskoordinate z, die die Nummer der jeweiligen Schicht 101, 102 und 103 wieder- 
gibt. In die Auswertung gehen die Werte des Abfragesignals S 1 und des Antwortsignals S2 an den zeitlichen Abtaststel- 
len t = k • X mit k = 1,2 . . . ein. 

[0043] Der eindimensionale Rekursionsalgorithmus basiert auf der tJberlegung, dass das eingespeiste Abfragesignal 45 
SI an jeder Grenzflache zwischen benachbarten Schichten 101, 102 und 103 teils reflektiert und teils transmittiert wird. 
Dadurch kommt es zu jedem diskreten Zeitpunkt und an jeder von dem Abfragesignal erfassten Grenzflache zu einer 
Aufteilung in einen vorlaufenden und einen zuriicklaufenden Signalanteil. Entsprechend der Nomenklatur der \Sfellenlei- 
tertheorie sind die vorlaufenden Anteile mit a und die zuriicklaufenden Anteile mit b bezeichnet. Beide Signalanteile a 
und b hangen von der diskreten Ortskoordinate z und von der diskreten Zeitkoordinate k ab. Vbrausgesetzt, dass die 50 
Sonde 1 und das iiberwachende Medium 30 zu Beginn der Messung energiefrei sind, ist das Antwortsignal S 1 zum ersten 
Betrachtungszeitpunkt k = 1 gleich dem am Sondenanfang reflektierten Anteil b(0,l) des in die Sonde 10 eingespeisten 
Abfrageimpulses a(0,0). Der Zusammenhang zwischen dem eingespeisten Signal a(0,0) und dem reflektierten Signal 
b(0,l) wird durch den Reflexionsfaktor R(l) an der Grenzflache zur ersten Schicht 101 mit der diskreten Ortskoordinate z 
- 1 bestimmt. 55 
[0044] Der Reflexionsfaktor R(z) an einer Grenzflache ergibt sich aus den Impedanzen Z(z) und Z(z-l) der beiden an- 
einandergrenzenden Schichten mit der diskreten Ortskoordinate z bzw. z-1: 



wobei sich die Impedanz Z(z) aus einem vorher berechenbaren konstanten Anteil Ag«,, der durch die Geometrie der ein- 
gesetzten Sonde bestimmt ist, sowie einem durch das zu iiberwachende Medium 30 bestimmten Anteil zusammensetzt: 
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[0045] Hierbei steht uO fur die natiirliche Permeabilitat, eO fur die natiirliche Dielektrizitatskonstante und er fiir die re- 
lative Dielektrizitatskonstante. Die relative Dielektrizitatskonstante er ist als einzige der eingehenden GroBen von der 
Ortskoordinate z abhangig, Gleichung (2) gilt insbesondere dann, wenn sich in dem Medium 30 nur dielektrische Kon- 
traste ausbilden. Das Medium 30 weist dann kein magnetisches Verhalten auf (d. h. u = uO). Dies ist bei praktisch alien 
gangigen Kuhlmitteln, d. h. auch bei Wasser, der Fall. Obergange zwischen verschiedenen Aggregatszustanden fuhren 
dann nur zu einer Veranderung der Dielektrizitatskonstante. Bei einem derartigen Medium 30 variieren die elektrische 
Impedanz Z(z) und eine Ausbreitungsgeschwindigkeit c(z) in der zugehorigen Schicht 101, 102 oder 103 nicht unabhan- 
gig voneinander. Es gilt: 

mo 



15 



20 



[0046] Damit lasst sich die Ausbreitungsgeschwindigkeit c(z) der aktuellen Schicht mit der Ortskoordinate z aus der 
Ausbreitungsgeschwindigkeit c(z-l) der vorhergehenden Schicht mit der Ortskoordinate z-1 sowie dem Reflexionsfak- 
tor R(z) der aktuellen Schicht berechnen: 



(4) 



25 



[0047] Ebenso kann auch die relative Dielektrizitatskonstante £r(z) der aktuellen Schicht aus der relativen Dielektrizi- 
tatskonstante er(z-l) der vorhergehenden Schicht sowie dem aktuellen Reflexionsfaktor R(z) der aktuellen Schicht be- 
rechnet werden: 



£r{z) = sr{z-\y — - ^ 



(5) 



30 [0048] Beriicksichtigt man auBerdem den Zusammenhang zwischen dem Reflexionsfaktor R(z) und dem Transmissi- 
onsfaktor T(z), der den eine Grenzschicht passierenden Anteil beschreibt: 

T(z) = J\-R(zf (6) 

35 sowie die Inversbeziehung zu Gleichung (1): 

*>-*-*£|!j 
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liegen alle Voraussetzungen vor, urn liber eine in Zeit und Raum sukzessive fortschreitende Betrachtung alle interessie- 
renden GroBen des Mediums 30 zu berechnen. Dies sind beispielsweise die Reflexionsfaktoren R(z), die Impedanz Z(z), 
die Ausbreitungsgeschwindigkeiten c(z), die relativen Dielektrizitatskonstanten er(z) aller Schichten 101, 102 und 103 
ebenso wie die ortlichen Positionen d(z) der Grenzflachen zwischen den benachbarten Schichten 101, 102 oder 103. Die 
45 Lage dieser Grenzflachen ergibt sich aus: 



d(z) = d(z-l)+T • c(z) (8) 

[0049] Das Prinzip wird anhand der in der Auswerteeinheit 3 fur die beiden ersten Abtastwerte des Antwortsignals S2 
50 durchgefuhrte Auswertung erlautert. Der erste Reflexionsfaktor R(l) wird aus den ersten Abtastwerten a(0,0) und b(0,l) 
des Abfragesignals SI bzw. des Antwortsignals S2 ermittelt. Unter Berucksichtigung der bekannten Impedanz der der 
Sonde 10 vorgeschalteten Vorlaufstrecke 15 sowie der Gleichungen (7) und (3) kann damit auch die Ausbreitungsge- 
schwindigkeit c(l) innerhalb der ersten Schicht 101 errechnet werden. Die Ausdehnung der ersten Schicht 101 und damit 
die Position der Grenzflache zwischen den Schichten 101 und 102 lassen sich dann uber Gleichung (8) bestimmen. 
55 Durch Zuordnung des ermittelten Reflexionsfaktors R(l) zu dem Ort des Sondeanfangs, d. h. des Anfangs der ersten 
Schicht 101, erhalt man den ersten Wert der gesuchten lokalen Verteilung des Reflexionsfaktorverlaufs R langs der 
Sonde 10. 

[0050] Der zweite Abtastwert b(0,2) des Abfragesignals S2 ergibt sich aus einer Summation des am Sondenanfang re- 
flektierten Anteils des zweiten Abtastwerts a(0,l) des Abfragesignals SI und des Anteiis des ersten Abtastwerts a(0,0) 
60 des Abfragesignals SI, der zunachst am Sondenanfang transmittiert, dann an der Grenzflache zwischen den Schichten 
101 und 102 reflektiert und am Sondenanfang wieder transmittiert worden ist. 

[0051] Wie auch aus der in Fig. 2 angegebenen Beziehung fur b(0,2) hervorgeht, ist die einzige unbekannte GroBe der 
Reflexionsfaktor R(2) an der Grenzflache zwischen den Schichten 101 und 102. Damit kann in dem zweiten Rekursions- 
schritt zunachst der Reflexionsfaktor R(2) ermittelt werden. Die weiteren interessierenden GroBen ergeben sich dann 
65 analog zu der fiir die erste Schicht 101 angegebenen Vorgehensweise unter Berucksichtigung der Gleichungen (2) bis (8). 
[0052] Bei den folgenden Rekursionsschritten werden jeweils alle neu dazukommenden vor- und zuriicklaufenden An- 
teile a(z, k) bzw. b(z, k) bestimmt, urn daraus und aus den bereits vorher bestimmten GroBen den Reflexionsfaktor R(z) 
und alle anderen interessierenden GroBen der aktuellen Schicht zu ermitteln. 
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[0053] Nach vollstandigem Durchlaufen des Rekursionsalgorithmus liegt dann der Verlauf des Reflexionsfaktors R(z) 
in Langsrichtung der Sonde 10 vor. In Fig. 1 ist schematised der auf diese Art und Weise ermittelte Reflexionsfaktorver- 
lauf R in Abhangigkeit der Ortskoordinate d ebenso gezeigt wie das anregende Abfragesignal SI und das empfangene 
Antwortsignal S2. Um aus den errechneten diskreten Werten des Reflexionsfaktors R(z) auf eine kontinuierliche Kurve 
zu gelangen, kann eines der bekannten Interpolationsverfahren verwendet werden. 5 
[0054] AuBer der in den Fig. 1 und 2 dargestellten Anregung durch ein Abfragesignal SI in Form einer geschalteten 
Gleichspannung, d. h. in Form einer Sprungfunktion, sind auch andere Anregungsfunktionen moglich. Bei spiels weise 
kann das Abfragesignal SI auch die Gestalt eines Chirp-Signals, eines Impuls-Signals oder eines Burst-Signals (= ge- 
schaltetes Sinus-Signal) annehmen. Die Sendeleistung des Abfragesignals SI liegt typischenveise im Bereich einiger 
mW. Die erzielbare Ortsauflosung wird im Wesentlichen durch die Bandbreite des Abfragesignals S 1 bestimmt. Deshalb 10 
wird iiblicherweise eine Bandbreite von mindestens 50 MHz verwendet. Typische Werte fur die eingesetzte Bandbreite 
liegen bei 1 GHz. Wenn es sich bei dem zu untersuchenden Medium 30 um Luft handelt, wird mit dieser Bandbreite eine 
Ortsauflosung von etwa 7,5 cm erreicht. Bei einem anderen Medium 30 mit einer hoheren relativen Dielektrizitatskon- 
stante er ist sogar eine noch bessere Ortsauflosung zu erreichen. Fiir Wasser oder Wasserdampf liegt dieser Wert unter 
1 cm. Weiterhin kann die Ortsauflosung beeinflusst werden durch die fur die Auswertung verwendete fiktive Laufzeit x. 15 
Solange die Abtasttheorie erfiillt ist und auch das Sendesignal eine ausreichende Bandbreite hat, lasst sich durch eine Re- 
duzierung der verwendeten fiktiven Laufzeit Er die Ortsauflosung und damit die Messgenauigkeit des Verfahrens erho- 
hen. 

[0055] GemaB Gleichung (5) kann aus dem Reflexionsfaktorverlauf R auch der Verlauf der relativen Dielektrizitats- 
konstante er in Langsrichtung der Sonde 10 ermittelt werden. Dies ist insbesondere auch deshalb von Interesse, da aus 20 
der Dielektrizitatskonstante er die Dichte p des Mediums 30 ermittelt werden kann. Fur Wasser und Wasserdampf als 
Medium 30 gibt es einen Zusammenhang zwischen der relativen Dielektrizitatskonstante er und der Dichte p sowie der 
Temperatur 0: 



er = l + ^-p+\^+A3+A4 S • p' +\'-^+A6-&+Al \- p* + 
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wobei Al bis A10 Konstante sind und # sowie p eine auf einen Bezugswert normierte Temperatur bzw. Dichte darstel- 
len. Gleichung (9) sowie die Konstanten A 1 bis A 10 sind aus dem Fachbuch "Properties of Water and Steam in Si-Units" 30 
von E. Schmidt und U. Grigull, Springer- Verlag, 3. Auflage, 1982, Tabelle 9, S. 175 und 176 bekannt. Fiir ein anderes 
Medium 30 lasst sich ein ahnlicher Zusammenhang zwischen der relativen Dielektrizitatskonstanten er und der Dichte p 
angeben. 

[0056] Gleichung (9) gilt allgemein. Bei einem eindeutigen Aggregatszustand hangen jedoch zusatzlich die Dichte p 
und die Temperatur # voneinander ab. Diese Abhangigkeit liegt sowohl bei fliissigem Wasser als auch bei Wasserdampf 35 
vor. Der Zusammenhang von Gleichung (9) ist in Fig. 3 fiir jeweils eindeutige Aggregatszustande dargestellt. Der gestri- 
chelt eingetragene Kurvenabschnitt320 ist der dampfformigen Phase und der mit durchgezogener Linie dargestellte Kur- 
venabschnitt 310 der flussigen Phase von Wasser zugeordnet. Die StoBstelle beider Kurvenabschnitte 310 und 320 be- 
zeichnet einen Zustand, bei dem sich der dampfformige und der fliissige Aggregatszustand annahemd gleich verhalten 
und auch nebeneinander in Form einer Mischphase vorliegen. Dies tritt bei der sogenannten kritischen Temperatur (& - 40 
374°Q auf. 

[0057] Anhand der eindeutigen Beziehung gemaB Gleichung (9) kann aus dem ermittelten Verlauf der relativen Di- 
elektrizitatskonstante er zum einen auf das Vorliegen eines dampfrormigen oder eines flussigen Aggregatszustands und 
zum anderen auf die Dichte p geschlossen werden. Hierbei ist wiederum die Ermittlung eines Dichteverlaufs in Langs- 
richtung der Sonde 10 moglich. Dies stellt verglichen mit einem Standard TDR- Verf ahren eine wesentliche Erweiterung 45 
der Messmoglichkeiten dar. 

[0058] Anhand des Dichteverlaufs p langs der Sonde 10 kann auch ohne weiteres eine Fiillstandsmenge oder eine Fiill- 
standshohe des Mediums 30 oder der Teilbefullung 31 im Behalter 20 ermittelt werden. Dies ist eine wichtige Diagnos- 
einformation fiir die Betriebssicherheit des Behalters 20. 

[0059] In Fig. 4 ist eine Sonde 11 dargestellt, die beispielsweise in der Vorrichtung gemaB Fig. 1 eingesetzt wird. Die 50 
Sonde 11 erstreckt sich praktisch iiber die komplette Lange des Behalters 20 also vom Behaiteroberteil bis zum Behal- 
temnterteil. Die Sonde 11 ist aufgebaut als koaxialer Wellenieiter mit einem AuBenleiter 110 sowie einem Innenleiter 
120. Der bei dem vorstehend beschriebenen Verfahren ermittelte Reflexionsfaktor R und die Impedanz Z wird bei dieser 
Ausgestaltung der Sonde 11 durch das Medium bestimmt, das sich im Zwischenraum zwischen dem Innen- und AuBen- 
leiter 120 bzw. 110 befindet. Damit das Medium 30 in diesen Innenraum eindringen kann, sind in dem AuBenleiter 110 55 
Offnungen 115 vorgesehen. 

[0060] Je nach Bedarf kann auch das untere Ende der Sonde 11 mit einer Offnung versehen sein. Der vorab einmal zu 
ermittelnde Geometrieparameter Ac c o von Gleichung (2) wird bei der koaxialen Geometrie der Sonde 11 im Wesendi- 
chen durch den AuBenleiterdurchmesser 111 sowie den Innenleiterdurchmesser 121 bestimmt. 

[0061] Andere Ausfuhrungsformen fiir die Sonde 10 sind ebenfalls moglich. Beispielsweise kann die Sonde auch als 60 
Eindrahdeiter vor der metallischen Behalterwand oder als Zweidrahtleiter oder als Hohlleiter ausgebildet sein. 
[0062] GemaB Gleichung (9) geht in die Ermittlung der Dichte p neben der relativen Dielektrizitatskonstanten er auch 
die Temperatur d ein. Zur Erhohung der Messgenauigkeit ist es deshalb von Nutzen, wenn neben dem Antwortsignal S2 
auch mindestens ein zusatzlicher Messwert fur die Temperatur im Behalter 20 erfasst wird. Die Genauigkeit der aus 
dem Abfragesignal S2 abgeleiteten MessgroBen wird um so hoher, je mehr und je bessere Informationen auch iiber die 65 
Temperaturverteilung langs der Sonde 10 vorliegen. Aus diesem Grund gibt es eine weitere Ausftihrungsform fiir eine 
Sonde 13, bei der mehrere Temperatursensoren 131 bis 135 an der Sonde 13 angeordnet sind. Diese weitere Ausftih- 
rungsform ist in Fig. 5 gezeigt. Sie enthalt temperaturabhangige Resonatoren aus einem insbesondere piezoelektrischcn 
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Material, wie z. B. Quarz. Fur eine Hochtemperaturanwendung mit Temperaturen von 800°C und hoher eignet sich ins- 
besondere Langasit als Material fiir die Temperatursensoren 131 bis 135. Die Sonde 13 enthalt auBerdem als weiteren 
Unterschied zur Sonde 11 von Fig. 4 einen hohlen Innenleiter 125. Beide Ausgestaltungen als Voll- oder als Hohlkorper 
eignen sich gleichermaBen fiir den Innenleiter 120 bzw. 125. 

5 [0063] Die Temperatursensoren 131 bis 135 konnen als Resonator oder als Laufzeitglied ausgebildet sein. Je nach Aus- 
fuhrungsform wird dann entweder eine temperaturabhangige Resonanzfrequenz oder eine temperaturabhangige Laufzeit 
fur die Temperatursensoren 131 bis 135 ermittelt. Beispielsweise erfolgt die Abfrage bei der Resonator-Ausfuhrungs- 
forrn iiber ein Sweep-Signal und bei der Laufzeitglied-Ausfuhrungsform im Impuls-Echoverfahren. Die Laufzeit inner- 
halb eines Laufzeitgliedes liegt dabei aufgrund der akustischen Ausbreitung deutiich uber der Laufzeit des Abfragesi- 

to gnals S 1 durch die komplette Sonde 13. 

[0064] Die Detektion des Antwortsignals S2 erfolgt in einer sehr kurzen Zeitspanne, vorzugsweise im Bereich zwi- 
schen einigen 100 ns und einigen us. Aufgrund dieser sehr kurzen Messdauer ist das beschriebene Verfahren auch un- 
empfindlich gegeniiber innerhalb des Mediums 30 auftretenden T\jrbulenzen. Die durch Tbrbulenzen hervorgerufenen 
Anderungen im Behalter und im Bereich um die Sonde 10 sind namlich stets deutiich langsamer als die fur die Messung 

15 benotigte Zeitspanne. Dabei spielt es auch Uberhaupt keine Rolle, dass die gegebenenfalls fur die Auswertung benotigte 
Zeit uber dieser kurzen Messdauer liegt, da das Antwortsignal S2 in der Auswerteeinheit 3 bei Bedarf zwischengespei- 
chert werden kann. 

[0065] Mitlels einer anderen Sonde 14 kann die Messgenauigkeit weiter erhdht werden. Diese spezielle Ausfuhrungs- 
form der Sonde 14 ist in Fig. 6 dargestellt. Die Sonde 14 enthalt in Langsrichtung mehrere Referenzmarken 141 bis 145, 
20 die bei Anregung durch das Abfragesignal S 1 zu einem Anteil S23 im Antwortsignal S2 fuhrt, der ebenfalls in Fig. 6 dar- 
gesteDt ist. Die Referenzmarken 141 bis 145 befinden sich an bekannten Stellen der Sonde 14. Bei korrekter Auswertung 
sollte dann ein Reflexionsfaktorverlauf R oder ein Impedanzverlauf Z diese Referenzmarken an den richtigen Orten wie- 
dergeben. Andernfalls konnen die erkannten Abweichungen zu einer Korrektur der Auswertung herangezogen werden. 

25 Patentanspriiche 

1. Verfahren zur Dichtebestimmung eines in einem Behalter (20) befindlichen Mediums (30) mittels einer in den 
Behalter (20) hineinragenden Sonde (10), bei dem 

a) ein Abfragesignal (SI) in die Sonde (10) eingespeist wird, 
30 b) ein durch das Abfragesignal (SI) hervorgerufenes Antwortsignal (S2) detektiert wird, 

c) ein Reflexionsfaktor- oder Impedanzverlauf (R, Z) des Mediums (30) langs der Sonde (10) mittels des Zeit- 
verlaufs des Antwortsignals (S2) und 

d) ein Dichteverlauf (p) des Mediums (30) langs der Sonde (10) mittels des Reflexionsfaktor- oder Impedanz- 
verlaufs (R, Z) bestimmt werden. 

35 2. Verfahren nach Anspruch 1, bei dem der Reflexionsfaktor- oder Impedanzverlauf (R, Z) rekursiv aus dem Ant- 

wortsignal (S2) ermittelt wird. 

3. Verfahren nach Anspruch 1 oder 2, bei zur Ermittlung des Reflexionsfaktor- oder Impedanzverlaufs (R, Z) 

a) das Medium (30) fiktiv in quer zur Sonde (10) verlaufende Schichten (101, 102, 103) mit jeweils gieicher 
vorgebbarer Laufzeit (T) des Abfragesignals (SI) durch die Schicht (101, 102, 103) eingeteilt wird und 

40 b) sukzessiv in Einspeiserichtung des Abfragesignals (SI) fur jede Schicht (101, 102, 103) 

bl) ein aktueller Reflexionsfaktor (R) an der Grenze zwischen einer aktuellen und der angrenzenden vor- 
hergehenden Schicht (101, 102, 103) mittels des Abfragesignals (SI), des Antwortsignals (S2) und der 
bereits ermittelten Reflexionsfaktoren (R) zwischen alien vorhergehenden Schichten (101, 102, 103), 
b2) eine aktuelle Ausbreitungsgeschwindigkeit (c) in der aktuellen Schicht (101, 102, 103) mittels der 

45 Ausbreitungsgeschwindigkeit (c) der angrenzenden vorhergehenden Schicht (101, 102, 103) und des ak- 

tuellen Reflexionsfaktors (R) und 

b3) eine in Einspeiserichtung des Abfragesignals (SI) gerichtete ortliche Ausdehnung der aktuellen 
Schicht (101, 102, 103) mittels der Laufzeit (%) und des aktuellen Reflexionsfaktors (R) bestimmt wird 
und 

50 b4) ein mittels der ortlichen Ausdehnung der angrenzenden vorhergehenden Schicht (101, 102, 103) be- 

stimmter Ort der Grenze zwischen aktueller und angrenzender vorhergehender Schicht (101, 102, 103) 
dem aktuellen Reflexionsfaktor (R) zugeordnet wird. 

4. Verfahren nach einem der vorhergehenden Anspriiche, bei dem zur Ermittlung des Dichteverlaufs (p) sukzessiv 
in Einspeiserichtung des Abfragesignals (SI) fur jede Schicht (101, 102, 103) 

55 a) eine aktuelle relative Dielektrizitatskonstante (er) mittels des aktuellen Reflexionsfaktors (R) und der rela- 

tiven Dielektrizitatskonstante (er) der angrenzenden vorhergehenden Schicht (101, 102, 103) und 

b) ein aktueller Dichtewert (p), der der aktuellen Schicht (101, 102, 103) zugeordnet wird, mittels der aktuel- 
len relativen Dielektrizitatskonstante (er) ermittelt wird. 

5. Verfahren nach einem der vorhergehenden Anspriiche, bei dem mittels des Reflexionsfaktor- oder Impedanzver- 
60 laufs (R, Z) ein Fullstand abgeleitet wird. 

6. Verfahren nach einem der vorhergehenden Anspriiche, bei dem mittels des Reflexionsfaktor- oder Impedanzver- 
laufs (R, Z) ein Aggregatszustand abgeleitet wird. 

7. Verfahren nach einem der vorhergehenden Anspriiche, bei dem als Abfragesignal (SI) ein sprungformiges Si- 
gnal, ein impulsfbrmiges Signal oder ein Chirp-Signal verwendet wird. 

65 8. Verfahren nach einem der vorhergehenden Anspriiche, bei dem das Abfragesignal (S 1) mit einer Bandbreite von 

mindestens 50 MHz, insbesondere von mindestens 100 MHz, eingespeist wird. 

9. Verfahren nach einem der vorhergehenden Anspriiche, bei dem die Temperatur innerhalb des Behalters (20) er- 
fasst wird. 
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10. Verfahren nach Anspruch 9, bei dem die Temperatur zumindest an einer S telle der Sonde (13) detektiert wind. 

11. Verfahren nach Anspruch 9 oder 10, bei dem die Temperatur zur Bestimraung der Dichte herangezogen wird. 

12. Verfahren nach einem der vorhergehenden Anspruche, bei dem iiber das Abfragesignal (SI) mindestens eine 
Referenzmarke (141, 142, 143, 144, 145), die an einer definierten S telle der Sonde (14) angebracht ist, abgefragt 
wird und ein durch die Referenzmarke (141 bis 145) im Antwortsignal verursachtes Teilsignal (S23) zur Ermitdung 5 
des Reflexionsfaktor- oder Impedanzverlaufs (R, Z) des Mediums (30) langs der Sonde (14) herangezogen wird. 

13. Vorrichtung zur Dichtebestimmung eines in einem Behalter (20) befindlichen Mediums (30) mittels einer in 
den Behalter (20) hineinragenden Sonde (10), umfassend mindestens 

a) Mittel (4, 5, 15) zur Erzeugung und Einspeisung eines Abfragesignals (SI) in die Sonde (10), 

b) Mittel (15, 3) zur Detektion und Auswertung eines durch das Abfragesignal (S 1 ) hervorgerufenen Antwort- 10 
signals (S2), wobei 

c) die Mittel (3) zur Detekdon und Auswertung ausgebildet sind zur Bestimmung eines Reflexionsfaktor- oder 
Impedanzverlaufs (R, Z) des Mediums (30) langs der Sonde (10) mittels des Zeitverlaufs des Antwortsignals 
(S2) und zur Bestimmung eines Dichteverlaufs (p) des Mediums (30) langs der Sonde (10) mittels des Refle- 
xionsfaktor- oder Impedanzverlaufs (R, Z). 15 

14. Vorrichtung nach Anspruch 1 3, bei der die Mittel (3) zur Detektion und Auswertung ausgebildet sind zur rekur- 
siven Bestimmung des Reflexionsfaktor- oder Impedanzverlaufs (R, Z). 

15. Vorrichtung nach Anspruch 13 oder 14, bei der die Mittel (4) zur Erzeugung ausgelegt sind zur Erzeugung eines 
Sprung-, impuls- oder chirpfbrrnigen Abfragesignals (SI). 

16. Vorrichtung nach einem der Anspruche 13 bis 15, bei der die Mittel (4) zur Erzeugung ausgelegt sind zur Er- 20 
zeugung eines Abfragesignals (SI) mit einer Bandbreite von mindestens 50 MHz, insbesondere von mindestens 
100 MHz. 

17. Vorrichtung nach einem der Anspruche 13 bis 16, bei der innerhalb des Behalters (20) mindestens ein Tempe- 
ratursensor (131-135) angeordnet ist. 

1 8. Vorrichtung nach Anspruch 17, bei der der mindestens eine Temperatursensor (131-135) an der Sonde (13) an- 25 
geordnet ist. 

19. Vorrichtung nach einem der Anspruche 13 bis 18, bei der mindestens eine Referenzmarke (141-145) an der 
Sonde (14) angeordnet ist. 
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